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We analyze topological charge contributions from classical SU(2) center vortices with shapes of 
planes and spheres using different topological charge definitions, namely the center vortex picture 
of topological charge, a discrete version of FF in the plaquette and hypercube definitions and 
the lattice index theorem. For the latter the zeromodes of the Dirac operator in the fundamental 
and adjoint representations using both the overlap and asqtad staggered fermion formulations are 
investigated. We find several problems for the individual definitions and discuss the discrepancies 
between the different topological charge definitions. Our results show that the interpretation of 
lattice quantities like FF in the background of center vortices is rather subtle. 



I. INTRODUCTION 

Since Savvidy [l[ we know that the QCD-vacuum is 
non-trivial and has magnetic properties. In lattice QCD 
it was shown 0] that vortices (quantized magnetic fluxes) 
condense in the vacuum and compress the electric flux 
between quark and antiquark to a string leading to con- 
finement. This center vortex model [S-Q seems to be 
a very promising candidate to explain the phenomena 
that dominate the infrared regime. Numerical simula- 
tions have indicated that vortices could also account for 
phenomena related to chiral symmetry, such as causing 
topological charge fluctuations and spontaneous chiral 



symmetry breaking (SCSB) [9h12||. In particular, [13 1 
states that the topological charge of a gauge field can be 
determined from the shape of P-vortices, i.e. from vortex 
intersections and writhing points. 

Center Vortices are based on a discrete gauge symme- 
try of the action. A non-trivial center transformation of 
all link variables in one time (or space) slice [2| 



U (x, t ) => zU (x, t ), z e Z 



N 



(1.1) 



leaves the action invariant. More generally, this transfor- 
mation can be expressed in terms of gauge transforma- 
tions on a periodic lattice 



U (x,t ) g(x,t)U (x,t )g r (x,t + 1) 



(1.2) 



which are periodic in the time direction only up to a Zn 
transformation: 



g(x,t + L t ) = zg(x,t ), z e Z 



N 



(1.3) 



This "singular" gauge transformation is not really a 
gauge transformation, since it affects of course the 
Polyakov loop, a gauge-invariant observable. Restrict- 
ing such a transformation to a finite volume of a slice, 
the three dimensional Dirac volume, increases the action 
by a surface contribution, the vortex action. However, 
center vortices also increase the entropy, compensating 
the rise of the action and give essential contributions to 



the vacuum configurations [14| . The appearance of link 
variables close to non-trivial center elements survives the 
continuum limit and leads to singular gauge fields. This 
implies the question, whether center vortices are lattice 
artifacts. One could give a positive answer, if removing 
such lattice artifacts would not influence QCD. But it is 
just the opposite, removing center vortices destroys con- 
finement and the topological charge vanishes 9]. This 
means center vortices are an essential ingredient of the 
QCD vacuum. 



II. TOPOLOGY ON THE LATTICE 

In lattice calculations there is a common method to 
determine the topological charge Q from the integral 
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16tH 



(2.1) 



where F^ v is expressed in terms of the plaquette field 

iV = U^{x)U v {x + fi)Ul(x + v)Ut{x). (2.2) 

This expression is derived in the continuum from the 
transition between vacua with different winding num- 
bers [IE El 



Q 



s 3 



(2.3) 



J M = ——=e^apTr{A v d a Ap + 2/3A u A a Ap). (2.4) 

57T 



Since F^F^ = 9 M J M , Q can be re-expressed as the 
above volume integral (|2.ip . On the lattice, continuity in 
space is lost and it seems that one should be able to view 
any lattice field configuration as being a discrete copy of 
a smooth continuum configuration. This would always 
be topologically trivial, since FF is a total derivative. 
Nonzero Q may come from field configurations contain- 
ing gauge singularities [l7l - [T9| . 
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Another possibility to analyze the topology of a gauge 
field is given by the Atiyah-Singer index theorem [20j- 
I23I ] . It states that the topological charge of a gauge field 
configuration is proportional to the index of the Dirac 
operator in thi s gauge field background. For the overlap 
Dirac operator [23l425l | in the fundamental representation 
the index is given by ind D[A] = n_ — ra+ = Q, where 
n_ and n+ are the number of left- and right-handed ze- 
romodes. The adjoint version of the index theorem reads 
ind D[A] = n_ - n + = 2NQ = 4Q, where N = 2 
is the number of colors and the additional factor 2 is 
due to the fact that the fermion is in a real represen- 
tation, hence the spectrum of the adjoint Dirac opera- 
tor iD is doubly degenerate. As described in [26| the 
improved staggered operator also produces eigenmodes 
which can clearly be identified as zeromodes and all re- 
sults in this paper show perfect agreement between the 
two fermion realizations, considering that the eigenvalues 
of the staggered fermion operator have a twofold degen- 
eracy due to a global charge conjugation symmetry in 
517(2). We therefore have ind D[A] = n_ - n + = 2Q 
for fundamental and ind D[A] = ri- — n + = 8Q for ad- 
joint (asqtad) staggered fermions. The lattice version of 
the index theorem is only valid as long as the gauge field 
is smooth enough and satisfies a so-called "admissibil- 
ity" condition. It requires that the plaquette values 
are bounded close to 1, the value for very smooth gauge 
fields. Sufficient, but not necessary bound for the "ad- 
missibility" of the gauge field are 1 1 — U^ IV \ \ < 1/30 [13] , 
or ||1 - U^W < [6(2 + 7(2))] = 0.04882 j^. 

In this paper we discuss a smooth lattice field con- 
figuration, which fulfills the admissibility condition and 
shows clearly a discrepancy between the integral of FF 
and the topological charge derived from the lattice index 
theorem. We work with thick, spherical vortices in SU(2) 
lattice gauge theory and extend our analysis of the topo- 
logical charge and the lattice index theorem. The prob- 
lem seems to be related to the singular nature of vortex 
configurations questioning the FF definition of topologi- 
cal charge. In fact it is incorporated to the SU(2) nature 
of our spherical vortex configuration, for which also the 
vortex picture of topological charge fails. Here we analyze 
this problem in more detail and start with the discussion 
of plane vortices in different Z7(l)-subgroups. 



III. PLANE CENTER VORTICES AND 
TOPOLOGICAL CHARGE CONTRIBUTIONS 

We investigated planar vortices parallel to two of the 
coordinate axes in SU(2) lattice gauge theory with re- 
spect to the location of Dirac zeromodes in [28]. The 
vortices are defined by links varying in a U(l) subgroup 
of SU(2), Uf_, — exp(i<^<Tj). The direction of the flux 
and the orientation of the vortices arc determined by 
the gradient of the angle <p, which we choose as a piece- 
wise linear function of the coordinate perpendicular to 
the vortex. Upon traversing a vortex sheet, the angle 4> 



increases or decreases by it within a finite thickness of the 
vortex. Center projection leads to a (thin) P-vortex at 
half the thickness [29| ■ If we consider these thick, planar 
vortices intersecting orthogonally, each intersection car- 
ries a topological charge with modulus |Q| = 1/2, whose 
sign depends on the relative orientation of the vortex 
fluxes [30j. The plaquette definition simply discretizes 
the continuum (Minkowski) expression of the Pontryagin 
index to a lattice (Euclidean) version of the topological 
charge definition: 

^ = 167r2 J d4xtx [^^^] 

= - ^~2 J d^xe^aptTlTapJ 7 ^] = J d A xE-B 

Our xy-vortices have only non-trivial zt-plaquettes, i.e. 
an electric field E z , while zt- vortices bear non-trivial xy- 
plaquettes corresponding to a magnetic field B z . The 
topological charge is then proportional to E Z B Z . If the 
angle for different vortex sheets rotates in the same 
U(l) subgroup then parallel crossings give Q = 1/2 and 
anti-parallel crossings give Q = —1/2. 

Constructing the vortex sheets out of different U(l) 
subgroups the individual crossings do not contribute to 
the topological charge if the color vectors of E- and B- 
plaquettes are orthogonal. Maximal center gauge still 
identifies two vortex pairs intersecting in four space-time 
points, but for parallel vortex pairs the topological charge 
from intersection points (Q=2) now differs from the FF 
definition. We find configurations with zero, one and 
even half topological charge. If <f> rotates e.g. in the o\- 
subgroup for the zy-vortex and in 02 for the zi-vortex 
then FF gives no contribution to the topological charge. 
Other configurations are shown in Fig. [TJ 

The action density in Fig. [TJl) already shows that such 
orthogonal color vector intersections might be suppressed 
due to higher action and in fact, at the intersection points 
we find maximally non-trivial (yt-) plaqucttcs 

P y t ~ {-io k ){-iai)i(j k iai 

These rough configurations also seem to trouble the Dirac 
operators, which do not find any zero modes. During 
cooling the Q = 1- and Q = 1/2-configurations are only 
meta-stable and soon turn to anti-parallel vortex pairs 
with vanishing topological charge, whereas the original 
Q = 0-configuration from two parallel vortex pairs in 
different J7(l)-subgroups prefers the 5 = 2Si ns t action 
minimum and ends up with parallel color vectors yield- 
ing Q — 2. For admissible gauge fields (after some cool- 
ing) we also find the correct numbers of zero modes and 
therefore the different definitions of topological charge 
agree with each other. Nevertheless, for the spherical 
vortex we immediately guarantee the admissibility con- 
dition but still get a discrepancy. 
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FIG. 1. Topological charge density for two parallel vortices 
intersecting in four points (same geometry as above). The 
first vortex sheet of the xy- vortex rotates <j> from zero to 7r in 
o i and on to 2ir in 02 for the second vortex sheet, whereas the 
zi-vortex only rotates in o\ (a) or starts with a 02 rotation 
and rotates the second vortex sheet in o\ (b) we get Q = 1 - 
or if the second zt- vortex sheet rotates in 0-3 we find Q = 1/2 
(below) distributed as shown in c) with an action density in 
the intersection plane shown in d). 
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FIG. 2. Topological charge during cooling the (meta-stable) 
Q = 1- and Q = 1/2-configurations. The intersecting vortex 
sheets end up with anti-parallel orientation and the topolog- 
ical charge vanishes. 



IV. THE "SPHERICAL VORTEX PROBLEM" 

The spherical vortex of radius R and thickness A was 
introduced in [3l| and analyzed in more detail in [26j . It 
is constructed with the following links: 



U (x ) = \ e *P ( [a (\r ~ r \)n ■ a) < = l,/i = 4 ^ 
11 elsewhere 



with n{r, t) 



r 

\r- f | 



(4.2) 
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FIG. 3. During cooling the Q — O-configuration the action 
reaches a plateau with S = 2Si n3 t and the parallel vortex pairs 
end up with parallel color vectors equivalent to the Q — 2- 
configuration. 



where r is the spatial part of x v and the profile function 
a is either one of aq.,a_, which are defined by 



a+(r) = 



'0 

f (I 



r < R 



r-R 



A 
2 



o_(r)= + 



R - f < r < R 
R+^<r 
r < R 

R-£ <r <R 







_ A 
2 

A 
2 

i?+#<r 



# , (4-3) 



4 • (4-4) 



This means that all links are equal to 1 except for the 
i-links in a single time-slice at fixed t = 1. The phase 
changes from to tt from inside to outside (or vice versa) . 
The graph of a_(r) is plotted in Fig. 2 in [3l|, giving 
a hedgehog-like configuration. In maximal center gauge 
and after center projection, this configuration shows a 
single, spherical vortex without any intersections and 
hence no topological charge. Since only links in the time 
direction are different from 1, the topological charge de- 
termined from any lattice version of FF vanishes for this 
spherical vortex configuration. The index of the con- 
sidered Dirac operators however is nonzero, resulting in 
Q = =pl, for a±. On a 136 3 x N t lattice the plaquettes for 
the spherical vortex, Eq. (|4.1[) . satisfy the "admissibility" 
condition. 

In order to understand the discrepancy we apply stan- 
dard cooling to the spherical vortex configuration. For 
many cooling steps, the index of the Dirac operator does 
not change, but the topological charge quickly rises close 
to =pl for a± while the action S reaches a (nonzero) 
plateau. So, the index of the overlap Dirac operator 
agrees with the topological charge after some cooling. In 
Fig. 2]we plot the cooling history for a spherical vortex on 
a 40 4 -lattice. For comparison we also plot the topological 
charge of an instanton during cooling, which looks pretty 
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FIG. 4. Cooling of a spherical vortex on a 40 4 lattice. The 
topological charge rises from zero to close to one for a = q_ 
(right scale) while the action S (in units of the one-instanton 
action Si ns t) reaches a plateau (left scale). For comparison we 
also plot the topological charge during cooling of an instanton. 

much the same as for the spherical vortex. In fact, the 
action and topological charge densities spread over more 
and more time-slices, developing a hyper-spherical dis- 
tribution like standard instantons. We conclude that our 
spherical vortex develops an instanton- like structure dur- 
ing cooling, in agreement with [32| , stating that the Haus- 
dorff dimension of regions where the topological charge 
is localized, changes towards the total space dimensions. 

However, the vortex structure of our initial configura- 
tion is removed after a few cooling steps, i.e. the spherical 
vortex shrinks very quickly. This is in agreement with the 
fact that a single instanton doesn't show any vortex con- 
tent [33j], but it clearly shows that cooling significantly 
changes contents of the initial gauge configuration. 

Further we apply some Monte Carlo steps to our spher- 
ical vortex configuration using the Metropolis algorithm 
with a small spread, i.e. adding small quantum fluctu- 
ations, and analyze the vortex structure and topological 
charge. Fig. [5] shows the action and topological charge 
during 100 Metropolis and another 100 cooling steps. 
The action rises during the Monte Carlo update and the 
spherical vortex percolates over the whole lattice while 
the topological charge fluctuates around zero. The index 
of the lattice Dirac operator however indicates topolog- 
ical charge Q = =pl for a±, and also cooling after the 
Monte Carlo update still reveals the correct result for Q 
via FF, showing the same behavior as in Fig. 2) 

We also want to emphasize the difference to the spher- 
ical vortex on an asymmetric lattice with time extent 
Nt = 2, which we analyzed in the second part of [26| . 
The vortex structure in one of the two time-slices sur- 
vives much longer during cooling and leads to a static, 
singular Dirac monopole before falling through the lat- 
tice. For completeness we also plot the cooling history of 
a spherical vortex on a 136 3 x 2 lattice in Fig- [HI The con- 
figuration initially satisfies the " admissibility" condition, 



FIG. 5. Monte Carlo (Metropolis) update and cooling of a 
spherical vortex on a 16 4 lattice. Even after 100 MC steps 
cooling reveals the correct topological charge Q via FF. 




o ' 1 1 1 1 1 1 o 

200 400 600 800 1000 1200 



cooling steps 

FIG. 6. Cooling of a spherical vortex on a 136 3 x 2 lattice. 
The topological charge definition FF initially gives Q — 0, 
in discrepancy with the index of the Dirac operator. Dur- 
ing cooling it rises to close to one for a = q_ resolving the 
discrepancy. The second plateau around 900 cooling steps in- 
dicates the fractional topological charge of a Dirac monopole, 
for more details see [2{J. 



but still shows the discrepancy between FF and the index 
of the Dirac operator. The second plateau around 900 
cooling steps indicates the fractional topological charge 
of the Dirac monopole, appearing only on the asymmet- 
ric lattice. On a symmetric lattice, we find a continuous 
transition to an instanton-like structure, and the vortex 
structure is lost. The singular gauge transformation is 
smoothed out in time direction at the cost of our cen- 
ter vortex. There seems to be no possibility of center 
vortices without gauge field singularities. This brings us 
back to the statement from section |TT] that topological 
charge on the lattice may need singularities of the gauge 
field, which seems to agree with the fact that removal of 
center vortices leads to vanishing topological charge. 
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V. DISCUSSION 

We finally want to resolve the problem of topologi- 
cal charge determination via FF for the above discussed 
spherical vortex configuration. The usual expression of 
topological charge Q (|2 . 1 1) does not take into account the 
periodic boundary conditions of the lattice, which mod- 
ify the dynamics considerably even in the large volume 
limit and lead to a different topological classification [§3] ■ 
The full expression for Q must contain possible twists in 
the boundary conditions allowed in the adjoint represen- 
tation [3514371 ]. However, such boundary twists may also 
be hidden in gauge singularities, especially in singular 
gauge transformations defining center vortices on the lat- 
tice. Usually one tries to evaluate the integrals for gauge 
invariant quantities like (I2.f I) in the axial gauge A = 0, 
where the singularities are transformed away [38]. How- 
ever, in our particular case of the spherical vortex there 
is no way to gauge-transform the singularity away, except 
by applying the inverse of the initial singular gauge trans- 
formation, which then would define the boundary twist, 
giving the only contribution to the topological charge of 
our trivial gauge field in A = gauge. 

An easy way of thinking how to determine the topolog- 
ical charge of a gauge field configuration is given in [1 8j j . 
Woit suggests to locate the gauge singularities, locally 
gauge transform them away using different gauges, and 
to measure the degrees of the maps relating the different 
gauges. The sum of these degrees will be a topological in- 
variant, the topological charge. For our configuration we 
easily may follow the instructions in the above reference 
by splitting the lattice into time slices. We now can only 
consider the time slice containing our spherical vortex, 
where the singular gauge transformation now defines the 
corresponding mapping and its degree the correct topo- 
logical charge. 

To state the solution of our problem in a mathematical 
manner, we consider the homotopy of the above spheri- 
cal vortex configuration, Eq. (|4.f [) . The i-links of these 
spherical vortices fix the holonomy of the time-like loops, 
defining a map Ut(x,t = 1) from the a;?/z-hyperplane at 
t = 1 to SU(2). Because of the periodic boundary condi- 
tions, the time-slice has the topology of a 3-torus. But, 
actually, we can identify all points in the exterior of the 3 
dimensional sphere since the links there are trivial. Thus 
the topology of the time-slice is R 3 U{oo} which is home- 
omorphic to S 3 . A map S 3 — > SU{2) is characterized by 
a winding number 

resulting in N = — f for positive and N = + 1 for negative 
spherical vortices. With this assignment the index of the 
Dirac operator and the topological charge after cooling 
coincide with this winding number N. Obviously such 
windings, given by the holonomy of the time-like loops of 
the spherical vortex, influence the index theorem (39l.l40|. 
which gives the correct definition of topological charge. 



Other, topologically motivated lattice constructions of 
Q from the gauge field are given in [TtJ and [4l|, where 
one compares the gauge rotations necessary in contiguous 
cells (hypercubes) to put each cell into the same (e.g. 
axial) gauge. This enables one to construct transition 
matrices v n<ll at the lattice sites n common to neighboring 
cells c(n) and c(n + fx) which can be used to derive a 
geometric definition of topological charge 

J p(n.fi.iy) 

+ / d 3 xTr[(v~* d v v n ^)(v~* d p v n ^)(v~* d„v n<fl )], 

where A denotes the lattice, p{n, \x, v) the plaquettes and 
/(n,/z) the faces (cubes) of a cell c(n). Evaluating this 
expression for our spherical vortex, the only non-trivial 
contribution is given by the second term for f(n,fi = 4) in 
the time-slice of our vortex, resulting in the expression for 
the winding number given above. All other terms vanish 
because of trivial transition functions or the vanishing of 
FF for our configuration. 



VI. CONCLUSIONS 

We reported on problems defining topological charge in 
the background of classical center vortices on the lattice. 
First, planar vortex sheets are constructed by U(l ) rota- 
tions in a way that they intersect orthogonally. These 
intersections are known to carry a topological charge 
Q = ±1/2, but if the gauge rotations are defined in dif- 
ferent Z7(l)-subgroups of SU(2), FF fails to detect the 
topological charge contributions. The use of the Dirac 
operator is not save in this case because of maximally 
non-trivial plaquettes, which seem to suppress such con- 
figurations in the functional integral. However, for "ad- 
missible" gauge configurations of classical, spherical cen- 
ter vortices we find a discrepancy between FF and the 
lattice index theorem, for both, overlap and asqtad stag- 
gered fermions in the fundamental and adjoint represen- 
tations. Numerically, the discrepancy equals the winding 
number of the spheres when they are regarded as maps 
S 3 -> SU{2). The problem arises due to the periodic 
boundary conditions on the lattice and the fact that cen- 
ter vortices are based on singular gauge transformations. 
The standard FF definition of topological charge is not 
applicable for such gauge fields. In our case we can re- 
gard the singular gauge transformation resulting in our 
spherical vortex as a boundary twist, leading to an extra 
contribution to FF, resolving the discrepancy. However, 
this result shows that the interpretation of topological 
charge via FF is rather subtle in the background of center 
vortices. The mentioned admissibility conditions do not 
guarantee that a naive plaquette or hypercube definition 
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of topological charge gives the correct result. This is only 
guaranteed in the continuum limit, where vortex configu- 
rations may become singular and validate the derivation 
of topological charge via FF. On lattice configurations 
the FF definition of topological charge should only be 
used after cooling, which however may change the gauge 
field content significantly. Besides geometrically moti- 
vated definitions of topological charge, the most reliable 
tool to study the topological structure of a lattice gauge 
field seems to be the Dirac operator. 
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